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Starting with a micropolar formulation, known to account for nonlocal microstructural effects 
at the continuum level, a generalized Langevin equation (GLE) for a particle, describing the pre¬ 
dominant motion of a localized region through a single displacement degree-of-freedom (DOF), is 
derived. The GLE features a memory dependent multiplicative or internal noise, which appears upon 
recognising that the micro-rotation variables possess randomness owing to an uncertainty principle. 

Unlike its classical version, the new GLE qualitatively reproduces the experimentally measured fluc¬ 
tuations in the steady-state mean square displacement of scattering centers in a polyvinyl alcohol 
slab. The origin of the fluctuations is traced to nonlocal spatial interactions within the continuum. 

A constraint equation, similar to a fluctuation dissipation theorem (FDT), is shown to statistically 
relate the internal noise to the other parameters in the GLE. 

PACS numbers: 46.05.-|-b, 05.40.-a, 05.40.Ca, 78.35.-l-c 


As the space/time length scale in the forcing or the 
triggered deformation mechanism becomes comparable 
with the internal length scale of the material, the system 
response at the macro-continuum scale is significantly in¬ 
fluenced by the material microstructure. This renders 
the classical continuum hypothesis of strictly local in¬ 
teractions, or its possible extension using adiabatic con¬ 
tinuation arguments, untenable in the modeling of such 
response. Non-local modeling techniques such as microp¬ 
olar [Ij, micromorphic [i| or gradient 0 theories, which 
aim at incorporating long-range inter-particle interac¬ 
tions, bring forth microstructural effects by introducing 
material length scales in the constitutive formulation. In 
materials like polymers, granular solids etc. where the 
length scale is of the macroscopic order, long-range ef¬ 
fects could predominate and, in such cases, predictions 
through nonlocal models, unlike the the classical contin¬ 
uum model, are in closer conformity with experimental 
observations. If a continuum model can be replaced by 
a collection of harmonic oscillators and the focus is on 
the motion of a small region represented by a particle 
(an oscillator) with a single predominant translational 
degree-of-freedom (DOF), one arrives at a generalized 
Langevin equation (GLE) for the DOF after including 
the coupling effects from the neighboring oscillators Q. 
The GLE, an expedient modeling tool that replaces the 
infinite dimensional continuum, is widely used in areas 
such as soft condensed matter physics and cell biology. 
Despite the correspondence as above between the GLE 
and the continuum, standard forms of the GLE do not 
include length scale information characteristic of nonlo¬ 
cal continuum theories and are therefore not equipped to 
describe the physically relevant microstructural effects. 

This work is partly motivated by the emerimental data 
shown in FigH] (adapted from Fig. 6 of [J]). These corre¬ 


spond to the mean square displacement (MSD) plots of 
temperature-induced Brownian particles in a polyvinyl 
alcohol (PVA) slab, extracted through light scattering 
studies: one corresponding to a local region marked by 
a focussed high-frequency ultrasound force and the other 
a global measurement valid also for the region. In the 
steady-state, the MSD in the ultrasound forcing case 
shows significant fluctuations, though over a much lower 
frequency band. These fluctuations are absent when the 
ultrasound forcing is switched off. Attempts at numerical 
simulations using the standard GLE under the known ul¬ 
trasound forcing fail to reproduce such fluctuations (see 
FigH]), the origin of which is believed to be in the long- 
range microstructural effects. 

In this letter, we aim to arrive at a GLE that carries 
such microstructural information. A key aspect is an in¬ 
ternal noise term in the GLE, which reflects an inherent 
randomness in the micro-rotations consequent to an un¬ 
certainty relation involving two strain operators of which 
only one contains micro-rotational information. Through 
the internal noise, evolution of the microstructural inter¬ 
actions, as manifested in the continuum scale, is charac¬ 
terized. Numerical simulations of the new GLE is shown 
to capture the steady-state fluctuations in the MSD. Ex¬ 
tending Kubo’s second fluctuation dissipation theorem 
(FDT) [51, a constraint equation, relating the damping 
memory kernel and intensities of both the initial thermal 
fluctuation and the internal noise, is derived. 

Discrete Hamiltonian from a micropolar perspective. — 
This work is founded on the premise that the adop¬ 
tion of micropolarity in the medium suffices to bring in 
the microstructural information. Accordingly, in pur¬ 
suance of the micropolar continuum theory, the de¬ 
formation kinematics requires each material point to 
have micro-rotation DOFs in addition to the standard 
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translational DOFs. The total mechanical energy func¬ 
tional (n) for a geometrically nonlinear micropolar body 
of material volume Vq may be represented as, If = 
5 Ivo + fvo Po'^dV. Here v is the 

linear velocity, u) the spin velocity, po the mass density, 
PqJI the micro-inertia tensor and ip the specific free en¬ 
ergy potential, all represented in Lagrangian coordinates. 
Assuming micro-rotations to be small, possibly an order 
smaller than the translation DOFs, a consistent spatial 
discretization of H, as elaborated in the supplementary 
material^ with the details of the symbols used, leads to 
the Hamiltonian (H) corresponding to a discrete repre¬ 
sentation of the body. In the discrete body, surrounded 
by a set of bath particles is a system particle whose pre¬ 
dominant translation DOF, also called the system DOF, 
is described by the OLE that we intend to arrive at. In¬ 
deed, except for the system DOF, all other DOFs appear¬ 
ing in H are henceforth referred to as the bath DOFs. 

Accordingly, we split the Hamiltonian as H — Hg + Hb, 
2 

where Hg = —h rcnnUn does not explicitly involve 

the bath DOFs and 
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Hamilton’s equations; details of the derivation are pro¬ 
vided in 0 . The equations of motion, in a compact form, 
are given by Eqn. m and (ED- 
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ue = {u’’" = {u 2 l ... USN dll ■■■ dsAfj'^consists 

of displacements and micro-rotations of the bath 
variables and Pue = {p 2 i ■■■ Psn Pen P93 n} , 
the corresponding linear and angular mo¬ 
menta. A is a diagonal matrix with nonzero 
entries { 1 /to 2 i ... l/msN l/(Ai) l/(/ 3 Ar)}- 
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. Ki, A and B are constant matri¬ 


ces of dimensions {6N — I) x (3A — 1), (3A — I) x 3N 
and 3A x 3Arespectively. g is a ( 6 A — I) dimensional 
constant vector. Omxn designates a zero matrix of 
dimension m x n and Wix( 6 Ar-i) in T is a constant 

matrix. Denoting {uj Pus}'^ as Yb, Eqn. ([2]) may be 
written as: 


Here i,j € {1,2,3} are the indices denoting the three 
cartesian coordinates. Summations over the repeated in¬ 
dices (i, j) are implied, rrika and Ika are respectively 
the mass and mass-moment of inertia of the a**' particle. 
The coordinate index fc, fictitiously introduced in m and 
I, serves to maintain indicial consistency with the vectors 
M, 0 and p appearing in the expressions and helps in us¬ 
ing Einstein’s summation convention in k. pia, peia are 
respectively the linear and angular momentum compo¬ 
nents and Uia, 0ia the displacement and micro-rotation 
components, all evaluated at the a*’*' particle in the 
direction. Of specific interest are un and pn, the dis¬ 
placement and the momentum components of the system 
particle in the required direction, for which the GLE will 
be written. For (j, i) = (I, 2), (2,3), (3,1), S = 1 and for 
other combinations of (j, f), S = — I. A is the number of 
particles in the discrete body. 

Formulation of the new GLE .— Using the discrete 
Hamiltonian FI, the governing dynamics for the system 
and bath variables, described in terms of the displace¬ 
ment DOFs and the momenta, are obtained through 


— KYf,-b ■iiii{0(gjv_i)xi S + Bd} (3) 


where 


g _ 0(6Af-l)x(6Af-l) 

[Ki K 2 ] 0(6Af-l)x(6Af-l). 

Multiplying Eqn. ([3]) with exp (—Kt) and integrating over 
[0,t], we get the following implicit expression: 


Yb (t) = exp (Kt) Yb (0) -b exp (Kt) . 


J exp (-Ks) Mil (s) {0(6Ar_i)xl g -b B 0 (s)}'^ds 


(4) 


Integration by parts on the term 
t _ X 

/exp (-Ks) Mil (s) {0(6Ar_i)xi g} ds of Eqn.dH) 
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would lead to an equivalent representation given in 
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Eqn.dll). 

Yb {t) = exp {Kt) Yb (0) + (t) {0(6Ar_i)xi g}^- 

exp (Ki) K~^uii (0) {0(6Ar_i)xi 9}^+ 

t 

exp(Kt) J exp (-Ks) K~^wii (s) {0(6Ar-i)xi 9 }^ds-\- 

0 

t 

exp(Kt) y exp (-Ks) uii (s) {0(6Ar-i)xi B6»(s)}"^fis 

0 

(5) 

Substituting Yb (t), as in Eqn.®, into Eqn.dD), we obtain 
the governing equations of motion (Eqn.®) for the sys¬ 
tem variables with the micro-rotation DOFs still tagged 
on. 

Niil 

dt \pii / 
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The constant element wq in Eqn.® additionally con¬ 
tributes to the stiffness due to micropolarity. This con¬ 
stant is the second element of the two-dimensional vector 
TK“^ {0(e7v_i)xi g] ■ 

Our interest is in deriving a GLE for the system 
DOF Mil, which requires eliminating the micro-rotation 
DOFs from Eqn.([S|). Inherent configurational uncer¬ 
tainty of the microstructure and its time evolution 
would seem to imply that the micro-rotation DOFs are 
treated as stochastic processes. A further justification 
of this viewpoint is provided through the noncommu¬ 
tativity of the symmetrized polar and nonpolar opera¬ 
tors, Ep 2 Y Ep ^ and Ejip ^ {^E^p -t- E^^ , 

wherein Ep = E and Enp = F'^F are respectively 
considered to be measures of polar and nonpolar strains. 
R and F are micro-rotation and deformation gradient 
tensors respectively @ . Specifically, we can arrive at the 
following Robertson-Schrodinger uncertainty Q relation 
involving these two operators. 




[A,B] is the commutator and {A,B} the anti¬ 
commutator of the operators A and B. For 
a given symmetric matrix operator C, (C) = 

/(X; X')^C(X')/(X; X')dX' defines the mean of 
C with respect to /. f{X;-) is a compactly sup¬ 

ported vector valued function of X' with an arbitrary 
support D containing X and must be normalized, i.e. 
II/IP = fdX' — 1. Thus f^f is interpretable 

as the density associated with a probability measure. 
Here, for instance, we choose / = ||^"^;|| , a normalized 
line segment. Such a mean, though a scalar quantity, 
should be interpreted as an operator wherever appropri¬ 
ate. a A = — (A))"'" (A — (A))^ is the variance (an 

uncertainty measure) associated with A. In terms of F 
and R, the inequality ([T]) may be recast as. 



(^R^F + F'^R^ (^F^F - ) 

( 8 ) 


The uncertainty relation is nontrivial only if the RHS of 
the inequality (jS)) is strictly greater than zero and this 
is indeed the case in general. The RHS of ([5]) can be 
zero only if PQf — 0 where P = R^F + F^R and 
Q = F^F — (^F^f'^. This result is obtained upon local¬ 
ization of the integral in the definition of (PQ) , given D 
is arbitrary. Since P is positive definite and Q generally 
nonsingular, PQ is nonsingular too and thus PQf yf 0. 
Clearly PQ is not skew symmetric and we may discount, 
almost surely, the other possibilities of PQf T f. This 
ensures, with probability 1, that the uncertainty rela¬ 
tion is nontrivial. In other words, while eliminating the 
micro-rotation DOFs, they should be treated as stochas¬ 
tic processes. 

Since, at a time instant, the mean of a micro-rotation 
DOF would typically be an order smaller than its trans¬ 
lational counterpart, we may approximately identify the 
micro-rotation DOF as a zero mean random variable. Re¬ 
taining such noise terms in the final GLE is then crucial, 
as nonlocality necessarily implies nondeterminism Q , an 
aspect generally overlooked in nonlocal continuum theo¬ 
ries. The randomness in the micro-rotations, consequent 
upon the uncertainty relation Eqn. 0, may be contrasted 
with that in the initial conditions due to thermal fluctu¬ 
ations, yielding an additive noise term in the GLE. This 
leads us to recast Eqn.® and write the GLE including 
the two noise sources as: 

fwiil 
dt \pii J 

J [g {s -1) uii (s) j U(i)J J 

0 0 

(9) 



where f) (•) is the memory kernel that may be found from 
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the second term on the right hand side of Eqn.®. ^ (t) 
is a linear combination of un (0) and elements of Yi, (0) 
at time t. An application of Lyapunov’s central limit 
theorem (CLT) yields ^ (t) to be a zero mean Gaussian 
random variable at time t. On similar lines, the internal 
(multiplicative) noise W*, arising from a weighted sum of 
the micro-rotations as seen from the last term of Eqn. (jG]) , 
may also be characterized as zero mean Gaussian for fixed 
s and t. We may rewrite the new GLE (including an 
external forcing term F (t) for completion) in its more 
familiar second order form. 


mil + ku+ / r] {s — t) ii (s) ds = 
Jo 


P (.t) + ^ (t) + 



Wl {s — t)u (s) ds 


( 10 ) 


In Ean. lfTOl) . m = mu denotes the mass, k = 
— -I- the stiffness and r] = —r} the damping 

memory kernel for the system variable. Note that, if mi- 
cropolar effects are not considered, i.e. and wq are 
identically zero, the usual form of GLE is retrieved. 

A fluctuation-dissipation (FD) constraint. — In addi¬ 
tion to the uncertainty constraint as in Eqn. ©, the 
causality condition may impose a constraint on the new 
GLE. The latter would be similar to the ED theorem typ¬ 
ically associated with the conventional GLE. We desig¬ 
nate the integral term on the right hand side of Eqn. m 
by y{t). A general scheme to represent W* could be 
via a Wiener chaos representation 3, 0. However, 
for illustrative purposes, we consider the special form 
W* = exp {a{s — t))C. (s) (a is some constant and ( (s) ^ 
M (0, cr (s))) that enables representing y (t) as the Marko¬ 
vian solution of a stochastic differential equation, a is the 
intensity of C- Thus Ean. dTUl) is equivalent to the coupled 
set of equations: 


mu-\-ku-\- / r] {s — t) ii (s) ds = F (t)f (t)y ft) (11) 
Jo 

y{t) =y{t)+u{t)C{t) (12) 

Under a random change of time t —> /3(t) [ni,M(t)C(t) 
may represented via a zero mean Brownian motion 

B {P (t)), where /3 (t) = ^cr (s)^ u (s)^'j ds. Eollowing 

the derivation of Kubo’s second EDT 3, the constraint 
equation may now be directly written as 

{{z (0) -b c (0)) (z (t) + f (t))) = T] (t) ksT (13) 
where z {t) = ^ exp {s — t) B (/3 (s)) ds^ . 


Some manipulations on Eqn. (USD would lead to 
Eqn. (IT4)) . 

J exp (s - t) (0) B {P (s))^ ds = r] ft) kBT-{^ (0) ^ (<)) 

(14) 

We see from Eqn. d that, as the left hand side vanishes 
in absence of micropolarity, we get back the usual FD 
theorem. 

Numerical simulations and experiments. — The numer¬ 
ical tests are carried out by simulating the standard GLE 
and the proposed one within a Monte Garlo setup. In 
both the simulations Prony series approximations are 
adopted for characterizing the memory kernels and the 
non-Markovian SDEs are converted to a set of Marko¬ 
vian SDEs [ 13 . From the nonlocal interaction effects, as 
anticipated, the fluctuations in the steady-state regime 
are indeed captured in the simulation via the proposed 
GLE (FiglU, which parallels the experimental observa¬ 
tion (FiglD). In FiglU we also see that simulation using 
the standard GLE fails to produce the steady-state fluc¬ 
tuations. The proposed GLE is further tested in the con¬ 
text of an inverse problem, wherein using it as a process 
model, a stochastic projection on the experimental MSD 
data through a nonlinear filter [l^ leads to an estimate 
quite close to the measurement (Figl3]). However, the 
same exercise with the standard GLE as process model 
(for the same Monte Garlo sample size) produces a com¬ 
pletely different response estimate. 



FIG. 1: Experimental plots of MSDs 

On a concluding note we emphasize that a suitably 
designed internal noise accounts for nonlocality, a fea¬ 
ture of space-separated multi-particle interaction, within 
a 1-dimensional evolution equation. The new GLE can 
more faithfully model the response of the material tak¬ 
ing into account some non-trivial aspects triggered by the 
microstructure. In future work, other than working out a 
more general model for the internal noise, it would be of 
interest to see if similar GLEs are derivable either from 
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2: Simulation via standard and proposed GLE 
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lattice models of materials or from Eshelby-type config¬ 
urational forces. 
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MICROPOLAR THEORY 

Micropolar theory [i[ , a generalized continuum approach, provides a phenomenological route to include length scale 
effects (generally manifested because of the micro-structural features of the material) in the response of a continuum 
to external stimuli. As opposed to the classical continuum model where the motion of each material point is described 
by its position vector alone, each such point in a micropolar continuum is additionally associated with a micro¬ 
structure which can rotate rigidly independent of its surrounding material. Therefore, every material point within the 
micropolar framework is endowed with six degrees of freedom (DOFs), namely three translational DOFs associated 
with the macro-element and three rotational DOFs with the micro-element. This is contrasted with the classical 
description wherein only the former three are present. Along with this difference in the kinematic description, the 
micropolar theory includes certain different features in the kinetic aspects too. The micropolar theory assumes the 
existence of a couple traction vector, in addition to the force traction considered in the classical continuum, for the force 
transfer mechanism across the interfaces of adjoining parts of the body. That is why, unlike the classical symmetric 
Cauchy stress tensor, a micropolar continuum generates an asymmetric force stress and a couple stress tensor. While 
the symmetric part of the force stress is responsible for the macro-deformation, the skew-symmetric portion along 
with the couple stress contribute to the rigid rotation of the micro-element. In the following three subsections, we 
briefly present the kinematics, equations of motion and constitutive relations for a micropolar continuum model. 


Kinematics 


Let Bq C be the reference configuration of a micropolar continuum body at time to- Macroscopically, each ma¬ 
terial point X of this continuum deforms to x in its spatial configuration, B, at time t. This translation is represented 
by a smooth map y as a; = y(X,t). y is a one-one and onto map; therefore possesses a unique inverse. As already 
mentioned, the microstructure associated with each material point of such a continuum can rotate independently 
along with the macroscopic deformation. To describe such rotational motion, assume an arbitrary vector H, at X in 
the reference configuration, which rotates to ^ at a; in its spatial configuration. In terms of a smooth invertible map 

X 

(, this motion is written as | = C (3j^)- The micro-rotation is described through a pseudo vector 0 = | 0 i 02 6 * 3 } 
with magnitude |0| := ip. Associated with the pseudo vector 6, there exists a proper orthogonal tensor R, using which 
the rotation map of the microstructure is written as ^ = R (0, t) H. A closed form expression for R in terms of 9 and 
p is given by the Euler- Rodrigues formula: 


R(6»,t) 


cos:/3l -I -W - 

P 


1 — cos p 

p^ 


e®e 


( 1 ) 


where I is the second order identity tensor, W the second order skew-symmetric tensor whose axial vector is 0 
(W = —e : 9, e being the third order alternate tensor or the Levi-Civita symbol) and 0 is the symbol for the exterior 
product or the dyadic product of vectors. 


Strain and Wryness Tensors 


Using the macro-translational motion y, the deformation gradient is defined as F = Similar to the well known 


polar decomposition in the classical theory, F could be uniquely decomposed as F = RU = VR. However a few 
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differences may be noted, viz. R here is not the macro-continuum rotation rather a measure of the micro-rotation, 
U and V, referred to as the right and left micropolar stretches respectively, are not symmetric unlike their classical 
counterparts. The right stretch tensor provides a material description of stretch whereas the left stretch tensor gives 
the spatial description of the same. For the definitions of the micropolar strain tensor E and the wryness tensor F in 
a Lagrangian description, the following may be adopted: E = U — I and F = — : (R'^VxR-)- The wryness tensor 

in the spatial configuration ( 7 ) is related to the material description as 7 = RFR"'’. 


Equations of Motion 


The physical requirement of balances of linear and angular momenta for the continuum must be met. Localization 
of these global balance requirements leads to the following equations of motion in the material description. 


Po ( ) = Vx ■ T -t- P 060 and 


( 2 ) 


j(0) =Vx-M + po^o-£:TFT (3) 

Here T and M respectively refer to the force stress and couple stress tensors. &o and /q are the externally applied body 
force and body couple densities respectively, po and J respectively refer to the mass density and the micro-inertia 
tensor, all expressed in the material coordinate. The tensorial form of the micro-inertia J may be simplified to po>^I 
in case of isotropic spin. Solution of the above equations of motion along with proper boundary and initial conditions 
would represent the response of the continuum. Owing to 0 evolving on the special orthogonal nonlinear manifold 
SO(3), integration of Eqn.([3]) is nontrivial. Specifically, the rotations of a material point at two different instances 
must be appropriately transformed to refer to the same tangent space of the manifold. However, for small 6, the 
necessary transformation is nearly the identity operator enabling integration within a purely linear space framework. 


Constitutive relations 


In the context of micropolar hyperelasticity, it is postulated that there exists a free energy functional (equivalently 
an internal strain energy functional considering only the mechanical response instead of the thermo-mechanical one) 
which defines the constitutive relations given as T = and M = Material description of the free 

energy functional is used to establish these relations. The specific form of the free energy density for the non-linear 
continuum considered here is given below. 

?/'(U,F,R) =^(U,R)-f V-(U,R)-k^(F,R) (4) 


where 

(U, R) = ip (tr (RUU'^R'^) - 3) -f (det (F)) 


poV- (U, R) = iff {tr (RUU'^R'^) - tr (RU^R"^)) 


(F,R) = i (RFR'^))^ -k ptr (RF^R^) + ^tr (RFF^rT)) 

The function h above represents an energy density associated with volumetric deformation, the functional form of 
which is /i = I ^det (F)^ — 1^ “ i In (det (F)) — pin (det (F)). Here p, A, it, d, /3 and 7 are the material parameters 
to be determined from experiments. 
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FORMULATION OF DISCRETE HAMILTONIAN 

The Hamiltonian (H) in a closed system may be given in terms of the total energy (H) comprised of the kinetic energy 
and the strain energy stored in the body. The kinetic energy of the body is given by K.E. = i fy^ po (v.v + Jlj.lj) dV 
and the potential energy stored is given by P.E. = fy^ pQipdV. These volume integrals are defined on the reference 

configuration where dV denotes a volume element and Vq the volume of the configuration Bq. is the 

linear velocity and a; = the spin velocity in the reference coordinate. Thus H = ^ Jy^ po {v.v + Jlj.uj) dV + 

fy^ PQipdV. As we consider micro-rotations and their gradients to be small, we may ignore the contribution of ip in 
E[. Moreover, the material considered (e.g. polymer) here is assumed to be elastic with nearly isochoric (det (F) 1) 

deformation. Hence, we may ignore the contribution of h too. Thus we arrive at a simpler description of H given as: 

H=l[ po{v.v + Jiv.iv)dV+f (lp{tr{FF'^) -3) + ]fi{tr{FF^)-tr{FIl^FIl^)))dV (5) 
2 Jvo Jvo V 2 4 J 

Towards writing H in terms of the displacement vector, u = y (X) — X, we use the relation, F = Vxu -f I and 
rewrite Eqn.dS]) as (for conciseness V is used instead of Vx)- 




Po {v.v + JiV.Lj) 



^-/itr ^Vu (Vu)'^-f Vu-f (Vu)'^^-I- ^ 

^atr ^Vu (Vu)"^ -f Vu -f (Vu)'^ ~ 

y^Ktr (VuR'^VuR'^ -I- VuR'^R'^ -k R'^VuR'^ -k R'^R'^)^ 


dV 


( 6 ) 


In what follows, we describe a consistent discretization of H given in Equ.®. Using the particle mass rriia = poa^Va 
and the micro-inertia lia = {poJ)^AVa for a = I,...,TV, the following discrete form of the kinetic energy may be 
arrived at. 


1 

2 


N 

{rriiavl 




(7) 


N is the total number of material particles. In mia and lia, the subscript i is introduced to maintain notational 
similarity with the other vector field variables (e.g. Via) and this is helpful in using Einstein’s summation convention 
in L In order to discretize the potential energy, a kernel approximation of the the displacement field is adopted. For 
example, we may write 

(tr(vu(Vu)'"))dF^^ [V x'<^> {X - X')]/.dV^dV, 

where <i> {y) is a sufficiently smooth symmetric kernel with respect to y. Uj is the j*’*' component of u . A 
subsequent particle representation leads to the following expression: 


X.(‘’'(^u(Vuk)) 


dV 


N 

E 

7=1 



N 


(y.,/3—1 


( 8 ) 


Here, Umr denotes the component of u evaluated at the material particle. 

= i:^({(Vx'4>)7}„{(V x^4^)Jj' AFaAVp^ AFj is the appropriate weight. Now, for small 0, we use the 
approximation R Ri I — e : 0 to get: 

(vu(RR)'^) dV= l^tr (R^VuRT) dV ^ Y + 2S (0^^)^^^-,)^^^. + [ 1 - 


tr 


I Vo 


2 

k'y 


Q!, 7 =l 


k^i 


( 9 ) 




4 


N 

where Wia-y = ({(^x'^’)7}q^^)^^7 - For {j,i) = (1,2), (2,3), (3,1), S = 1 and for other combinations of 

7=1 

S = — 1. Neglecting contributions from tr (VuR^VuRT) and tr(RR)'^, the discrete Hamiltonian finally takes the 
following form. 



\lTlict ^ioL J a t 1 

a — 1 ^ ' a,p — l a,7=1 




i^j 


( 10 ) 


EQUATIONS OF MOTION OF THE DISCRETE SYSTEM 


The governing dynamics for the system and bath variables, described through Hamilton’s equations in terms of the 
displacement DOFs and the momenta, are given by the equation pairs (nmni) and (dam respectively. 


— ( 1 - ZlL 

dt dpii mil 


( 11 ) 



dH 

duii 


N N 

WllllUii + ^ WiiipUi/S + ^ Sill/? 

/?=2 7=1 



( 12 ) 


d dH 

dt dpuB 


P21 

m2i 


Pen 

hi 



(13) 


d „ dH 

~jtPu9 — 
at OUQ 


N N 

Wil2jiUip + Wi2^Lgi2^ 

/?=! 7=1 

N N 

Wi3N0Uii3 + WiNjLesN^ 

/?=! 7=1 

Ce {uii...U3n)'^ + D {uii...U3n)'^^ 


(14) 


where Lei 2 -y = 


(di-fdi'j + 2S (^fc7);s/j,i) + (^fc2)i=i 


and Lg 3 N~f = 


(ds-ydij + 2S + {^kN)i= 


ug = = {u 2 i ... U 3 N dll dsAT j'^consists of displacements and micro-rotations of bath variables and 

PuB = {P 21 P 3 N PBii Pb 3 n}'^, the corresponding linear and angular momenta. Cg is a matrix each of whose 

elements is a linear combination of 9ia where i G {1,2, 3} and a G {1, 2,..., JV}. D is a matrix of constant coefficients. 

Siiii and Wiii /3 are constants related respectively to Wmi and Wmis through the constant multipliers resulting from 
the partial derivatives of the Hamintonian. A similar definition holds for the constants Wia-y. 
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